Astronomy & Astrophysics manuscript no. ms 
May 25, 2010 



©ESO2010 



Light-element Abundance Variations in the Milky Way Halo 



Sarah L. Martell and Eva K. Grebel 



o 



(N 

< 

o 

Oh 
6 



> 

O 

o 

O 
O 



X 



Astronomisches Rechen-Institut 

Zentrum fur Astronomie der Universitat Heidelberg 

69120 Heidelberg, Germany 

e-mail: martell@ari .uni -heidelberg.de 

Received 26 January 2010 / Accepted 18 May 2010 

ABSTRACT 

We present evidence for the contribution of high-mass globular clusters to the stellar halo of the Galaxy. Using SDSS-II/SEGUE 
spectra of over 1900 G- and K-type halo giants, we identify for the first time a subset of stars with CN bandstrengths significantly 
larger, and CH bandstrengths lower, than the majority of halo field stars, at fixed temperature and metallicity. Since CN bandstrength 
inhomogeneity and the usual attendant abundance variations are presently understood as a result of star formation in globular clusters, 
we interpret this subset of halo giants as a result of globular cluster dissolution into the Galactic halo. We find that 2.5% of our sample 
is CN-strong, and can infer based on recent models of globular cluster evolution that the fraction of halo field stars initially formed 
within globular clusters may be as large as 50%. 
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1. Introduction 

Hierarchical structure formation is presently the dominant ex- 
planation for galaxy formation, based on observed fluctuations 
in the cosmic microwave background (Hinshaw et al. 2007) and 
sophisticated numerical simulations of their evolution to the 
present day (e.g., Diemand et al. 2007; Springel et al. 2008). 
In this picture, galaxies like the Milky Way are formed through 
the coalescence of multiple low-mass galaxies which develop 
within a much larger dark matter halo. The initial disagreement 
between the calculated mass function of dark matter subhaloes in 
the simulations and the observed mass function of nearby dwarf 
galaxies (e.g., Klypin et al. 1999; Moore et al. 1999) is being ad- 
dressed from several directions at the same time, from complex 
semianalytic simulations that include star formation and feed- 
back processes, and calculate the chemodynamical evolution of 
Milky Way-like galaxies (e.g., Johnston et al. 2008; Tumlinson 
2010 and references therein) to searches for extremely low-mass 
galaxies in the Local Group (e.g., Zucker et al. 2006a, 2006b; 
Belokurov et al. 2007). 

The stellar halo of the Milky Way is thought to have been 
constructed mostly through the early (^ 10 Gyr ago) accre- 
tion of low-mass protogalaxies. The halo exhibits considerable 
substructure in density and in velocity (e.g. Bell et al. 2008), 
and kinematically distinct streams presently observed in the 
halo (e.g., Majewski et al. 2003; Duffau et al. 2006; Martinez- 
Delgado et al. 2007) are interpreted as remnants of more re- 
cent or on going merger activity. The "ECHOS" identified in 
Schlaufm an et alj (12009) are interpreted as older substructure 
that has lost some spatial coherence with time. Studies of the 
abundance distributions and star formation histories in nearby 
dwarf galaxies (e.g., Koch et al. 2007a, 2007b, 2008a, 2008b; 
Kirby et al. 2008; Aoki et al. 2009; Frebel et al. 2009, 2010) 
have shown that there is a reasonable concordance between the 
properties of the present-day Milky Way halo (as characterized 
by, e.g., Schoerck et al. 2009) and the dwarf galaxies that would 



have been available as stellar contributors early in Galactic his- 
tory (e.g., Font et al. 2006; Carollo et al. 2007). 

However, the ongoing dissolution of globular clusters such 
as Palomar 5 (e.g., Odenkirchen et al. 2001, 2002, 2003; 
Rock osi et al. 2002; Grillmair & Dionatos 2006) and NGC 
5466 (Belokur ov et alJ l2006) implies that some fraction of halo 
stars are initially formed in globular clusters. This dissolution is 
driven by internal 2-body relaxation, stellar evolution processes 
and tidal interactions with the Galactic potential, and can cause 
significant mass loss over the lifetime of typical halo globular 
clusters (e.g., Gnedin & Ostriker 1997). 

One particular mode l of globular cluster formation, de- 
scribed in D'E rcole et alJ (1200 8). posits that all old halo clusters 
surviving to the present day have lost at least 90% of their initial 
mass. Early in the development of the cluster, winds from AGB 
stars collected in the cluster center and formed a second gen- 
eration of low-mass stars. Type la supernovae then caused the 
cluster to expand, and stars at large radii, mostly members of the 
first generation, were lost. In the model, the end result, ^10 Gyr 
later, is a 10 5 - 10 6 Mq cluster with two stellar populations that 
differ slightly in age and abundance pattern. 

This two-generation model was developed specifically to ex- 
plain the anomalous light-element abundance patterns observed 
in globular cluster stars, specifically the presence in every old 
globular cluster in the Milky Way of a subpopulation with typ- 
ical Population II abundances, and a second subgroup with the 
same metallicity but enhanced N, Na, and Mg along with de- 
pleted C, O, and Al. This abundance bimodality has been stud- 
ied extensively in globular clusters (Langer et al. 1992, Kraft 
1994, and Gratton et al. 2004 all provide thorough reviews 
of the topic), and explanations fo r the second a bundance sub- 
group have varied from pollution (Cannon et al. 1998) to inter- 
nal mixing (Langer 1985) to enrichment of star-forming gas by 
moderate- to high-mass stars (e.g., Cottr ell & DaCosta 1981 ; 
Yong et al. 2008) or high-mass binaries ddeMinket al.ll2009l) . 
Surface pollution of already-formed stars would result in abun- 
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dance anomalies that are erased at first dredge-up, while current 
models of deep mixing (e.g., Charbonnel & Zahn 2007) indi- 
cate that it only begins to operate after first dredge-up. Both 
of these processes have been effectively ruled out as explana- 
tions for abundance bimodality by the presence of abundance 
variations at all evolutionary phases in globular clusters (e.g., 
Briley et al. 2002; Harbeck et al. 2003). In the primordial en- 
richment scenario, there is ongoing discussion over the exact 
source of enriching material. In some models the source is mod- 
erately high-ma ss (~ 4 — 5M(7)) AGB stars (e.g., Parmentier et 
al. 1999), while iDecressin et a l. (2007) claim that rotating high- 
mass (M S IOMq) stars are a better source for processed ma- 
terial because of their very short lifetimes and de Mink et all 
(2009) prefer high-mass binaries because of their potentially 
strong mass loss and low wind velocities. 

Stars with these light-element abundance anomalies are 
readily identified through strong UV/blue CN molecular absorp- 
tion and relatively weak absorption in the CH G band, and are 
hereafter called "CN-strong stars", with the understanding that 
the full abundance pattern from C through Al necessarily follows 
the CN variation. They are not observed to exist in open clus- 
ters (e.g., Smith & Norris 19 84; Jacobson et al. 2 008; Martell & 
Smith 2009) or the halo field dGratton et alj2004l) . This feedback 
process apparently only occurs in the high-density environment 
of globular clusters, and as a result the characteristic sawtooth 
abundance pattern from carbon through aluminium can be used 
as a marker of globular cluster origin. 

Given the contributions globular clusters are presently mak- 
ing to the halo field, and the significant mass loss predicted theo- 
retically over the lifetime of the Galactic globular cluster system 
(e.g., Baumgardt et al. 2008), it is intriguing that no CN-strong 
stars have to date been observed in the halo. We interpret this as 
a qualitative sign that the contributions to the halo field of glob- 
ular clusters as we know the m today are relatively minor, as is 
also suggested in lYong et alj d2008l) . 

To test this interpretation, we searched for CN-strong halo 
giants in the Sloan Extension for Galactic Und erstanding and 
Exploration (SEGUE) survey dYannv et alj|2009t) . The SEGUE 
survey is a spectroscopic ex tension of imagi ng taken during the 
Sloan Digital Sky Survey (Yorketal. 2000), with targets se- 
lected to address questions of halo substructure and Galactic 
formation history. Data were taken from 2005 August through 
2008 July using a 640-fiber multiobject spectrograph and the 
same telescope at Apache Point Observatory that was used for 
SDSS imaging. The first portion of the SEGUE data was made 
publically available in 2008 as part of SDSS Data Release 7 
(DR7), including roughly 240,000 spectra in 200 "pencil beam" 
lines of sight containing stars chosen for specific purposes (i.e., 
fields in the Sagittarius stream, M dwarfs to study extremely lo- 
cal kinematic substructure, G and K giants to study the distant 
halo). In addition to flux-calibrated spectra, DR7 also offers the 
products of the "SEGUE Stellar Parameters Pipeline" (SSPP), 
which include derived stellar parameters like effective temper- 
ature, [Fe/H] metallicity, and radial velocity, determined auto- 
matically through template matching, x 2 minimization, cross- 
correlation or grid-matc hing methods, as appropriate. Lee et al. 
(2008a; 2008b) and lAllende Prieto et alJ d2008h give thorough 
explanations of the SSPP pipeline and process. 



2. The Data Set 

Data were obtained from SDSS DR7, through the online Catalog 
Archive Server. We selected all SEGUE plates, and from those 
all stars with [Fe/H]< -1.0, log(g) < 4.0, (g - r) > 0.2, and 




Fig. 1. Hertzsprung-Russell diagram for the initial SEGUE data 
set, split into two subgroups. 5066 G and K giants are shown 
as small points, while density contours for the remaining 17718 
stars rejected based on (T e g, log(g), [Fe/H]) position are shown 
as solid lines. 



mean signal-to-noise per pixel (SNR) larger than 20 were cho- 
sen. Additionally, the errors on various derived parameters were 
required to be small: cri og ( g ) < 0.5, 0"[ Fe/ / H ] < 0.5 with at least 
three independent [Fe/H] determinations, and reduced^ 2 of the 
best-fit template spectrum less than 2.0 both in the region of the 
Ca II H and K lines and in the CH G band. 

This intentionally generous set of selection criteria then had 
to be further sub-selected to isolate halo giants. To accomplish 
this sub-selection, we lowered the limit on log(g) to 3.0, then 
divided the initial data set into 0.2-dex-wide bins in [Fe/H], and 
removed AGB, main-sequence and turnoff stars from the H-R 
diagram of each of those subsets by rejecting all points more 
than 3<x in (g - r)o from the fiducial sequence. We also made the 
SNR requirement more stringent, requiring that the mean SNR 
per pixel in the wavelength range 4000 < A < 4100 be larger 
than 15. This left 5066 halo giants, out of the original 22,784 
stars. In Fig.Q] the small points represent those halo giants, while 
density contours of stars rejected based on CMD position are 
shown as solid lines. 

To facilitate later analysis, we convert the dereddened appar- 
ent ro magnitudes given in the SEGUE data to absolute M, mag- 
nitudes through a simple photometric paral lax calculation. We 
created a grid of 10 Gyr Padova isochrones (Marigo et al. 2008) 
with metallicities ranging from [Fe/H]= —1.0 to [Fe/H] = -1.8 
(the lowest metallicity available) at a spacing of 0.1 dex, and 
interpolate between those to match the metallicity of each in- 
dividual star. The absolute M r magnitude corresponding to the 
observed (g - r)o color of the star, on the interpolated isochrone, 
is then assigned as the true absolute M r magnitude. 

Errors in M, were calculated by Monte Carlo sampling of un- 
correlated errors in SSPP metallicity, go magnitude and r^ mag- 
nitude, drawn randomly from a Gaussian with a width equal to 
the reported errors on those quantities. This addition of error was 
done 10 4 times for each star, and we take the standard deviation 
in those 10 4 determinations of M r as the error on M T . This er- 
ror has a typical value of 0.3 magnitudes, with the largest values 
(k 0.5 magnitudes) on the 3% of stars with the largest errors in 
apparent ro and go- The age of the isochrones used had minimal 
effects on the derived M r values, with a change of only +0.06 
magnitudes for a shift of +1 Gyr. The limited metallicity range 
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Fig. 2. CN bandstrength index 5 (3839) versus dereddened (g—r) 
color for the G and K giants indicated in Fig.Q] As expected, the 
stars are mainly CN-weak, and there is a slight temperature de- 
pendence. As mentioned in the text, the stars with dramatically 
high 5 (3839) are mainly CEMP and carbon stars (shown as open 
triangles). 



of the isochrones is sufficient for the purposes of the analysis 
in Section 3, but lower-metallicity isochrones would allow us 
to study the distance distribution of the roughly 1/3 of our final 
dataset at lower metallicity, to compare the spatial distribution 
of our final data se t to the inner and outer halos identified in 
ICarolloetalJd200l . 

We measured 5(3839) dNorriset alJll98lh . a bandstrength 
index for the CN band at 3883A, for all halo giant spectra. 
5(3839) measures the magnitude difference between the inte- 
grated flux in the CN feature and the integrated flux in a nearby 
continuum band, with more absorption in the feature resulting in 
larger bandstrength. As can be seen in Fig. [2 there is a strong 
concentration at low 5(3839) in our data set. This is to be ex- 
pected, since the halo is primarily composed of CN-weak stars 
with typical Pop. II abundances. There is also a clear trend with 
temperature, in the sense that cooler stars have larger CN band- 
strengths. The cooling of stars as they ascend the giant branch 
reddens the spectra and permits more CN molecule formation, 
both of which increase the flux difference between the feature 
and continuum bands of 5(3839). There are, however, interest- 
ing outliers in Fig. [2J many of the stars with dramatically large 
5(3839) are carbon stars (shown as open triangles), with cor- 
respondingly large CH and C2 bandstrengths. Figure [3] shows 
four sample spectra from the halo giant data set: the uppermost 
spectrum (of SDSS J035 123.90+09245 1.3) is a low-metallicity 
carbon star, the next (of SDSS J 1 15934.87+002748.0) is a possi- 
ble CEMP star (Carbon-Enhanced Metal-Poor, having [Fe/H]£ 
-2.0 and [C/Fe]> +1.0, and described in Lucatello et al. 2006 
and references therein), with strong CH absorption redward of 
the G band and a low metallicity, the next spectrum (of SDSS 
J06441 1.96+275351.6) is not a carbon star, but is CN-strong, 
and the lowest spectrum (of SDSS J145301.24-001954.1) is a 
typical CN-weak halo star. 

The broad molecular features in the example carbon-star 
spectrum are quite clear. We use the strength of the CH feature 
around 4350A and the Swan (1,0) C 2 band at 4737A to identify 
carbon and CEMP stars and remove them from the final data set. 
Specifically, we measure these indices: 



Fig. 3. Sample spectra from the four main types of star com- 
prising our G and K giant data set: the topmost is from a car- 
bon star, with clear CH and C2 bands, and the next is a possible 
CEMP star, with a metallicity of [Fe/H]= -2.2 and strong CH 
absorption. The next spectrum down is from a CN-strong star, 
with stronger absorption in the 3883A CN band than the low- 
est spectrum, which shows the weak UV CN absorption seen in 
98% of our halo giant sample. The vertical lines mark the edges 
of the 5 (3839) bands: feature (3846A to 3883 A) and continuum 
(3883Ato3916A). 
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and consider all stars with [Fe/H]< -1.8 and s(c0) > -0.093; 
-1.8 <[Fe/H]< -1.4, s(c0) > -0.05, and s(cl) > 0.15; and 
[Fe/H]> -1.4, s{cQ) > -0.02, and s(cl) > 0.158 to have "strong 
carbon features". Although these carbon and CEMP stars tend to 
have strong UV CN bands, as is demonstrated in, e.g jAoki et al.l 
(120021) . that is a result of the unusually large ratio of carbon to 
oxygen in their atmospheres, and is not a result of the anticor- 
related C-N abundance pattern observed in CN-strong globular 
cluster stars. There are 109 stars removed from the data set be- 
cause of strong carbon features, but CN-strong stars are not en- 
tirely eliminated: carbon and CEMP stars are shown as open tri- 
angles in Fig. [2] and there are clearly other stars with stronger 
5(3839) than the main group, at fixed (g - r)o color. Of the 
109 stars we identify as carbon or CEMP stars, 89 have metal- 
licities below [Fe/H]= -1.8, 11 have metallicities in the range 
-1.8 <[Fe/H]< -1.4, and 9 have a metallicity above [Fe/H] = 
-1.4. This is not a s urprising metallicity d istribution: many stud- 
ies, including, e.g., Groenewegen ( 19991). iLucatello et aD (12006) 
and lAoki et al.l d2008l) . have found that the proportion of carbon 
stars rises at low metallicity. 



3. Distribution of CN and CH bandstrengths 

Since the purpose of this study is identifying stars in the halo 
with the CN-strong chemical signature associated with globular 
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Fig. 4. 5(3839) versus M T and a generalized histogram of 
55(3839) for M3 red giants observed by SEGUE. The classic 
bimodal CN distribution is apparent. 



clusters, and our data source is the moderate-resolution spec- 
troscopy provided by SEGUE, our methods follow closely the 
techniques used in low-resolution spectroscopic studies of CN 
bandstrength behavior in globular clusters. For example, the 
spectral feature used is the same 3883A CN band, and the index 
used to measure bandstrength is the well-known 5(3839). One 
first check of this approach is comparing CN behavior in well- 
studied globular clusters to halo stars with similar metallicities 
observed as part of SEGUE. Figure |4] shows two views of the 
CN bandstrength data for M3, measured from SEGUE spectra of 
M3 giants (membership information from J. Smolinski, private 
communication). The left panel shows 5(3839) versus absolute 
M r magnitude, and the typical globular cluster pattern is read- 
ily visible: the 5(3839) distribution separates into two parallel 
loci, rising with increasing luminosity. CN-strong stars (shown 
as open circles) are clearly distinct from CN-weak stars (filled 
circles) at fixed luminosity. 

The right panel of Fig. |4]converts the raw 5 (3839) data into 
a generalized histogram. In order to remove the temperature- 
related trend visible in the left panel, we fit a line to the CN- 
weak locus and measure the quantity 55(3839), the vertical dif- 
ference in 5(3839) between every point and the baseline. We 
then draw a generalized histogram by representing each point 
as a Gaussian centered at 55(3839), with a FWHM equal to the 
error on 55(3839), and summing the individual Gaussians to- 
gether. In Fig. |U we base the Gaussian widths in the general- 
ized histogram on the actual measurement errors on 5 (3839), as 
determined by Monte Carlo sampling of the spectral error vec- 
tors reported by SEGUE. Typical values of crS (3839) are 0.014 
magnitudes, and we amplify that error by a factor of four in the 
construction of the generalized histogram. This allows for errors 
not accounted for in our Monte Carlo noise sampling, particu- 
larly errors in flux calibration and the SSPP-derived parameters. 
We calculate generalized histograms for the CN-strong and CN- 
weak stars independently, and normalize both curves to the peak 
of the CN-weak curve, so that the relative heights represent the 
relative numbers of stars in the two groups. 

Figure[5]shows analogous data for field stars in the final data 
set with [Fe/H] within +0.05 dex of M3, and similar structure 
can be seen to Fig. |4] with interesting differences in scale. As in 
Fig. [4] we fit a baseline to the CN-weak stars in the M, - -5 (3839) 
plane (left panel), measure the quantity 55(3839), and convert 
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Fig. 5. 5(3839) versus M r and a generalized histogram of 
55(3839) for field red giants with metallicities similar to M3 
([Fe/H] between -1.62 and -1.52). Although CN-weak stars 
comprise a much larger fraction of the overall data set, there are 
still relatively CN-strong stars seen. The CN-strong curve has 
been magnified by a factor of 30 for clarity. 



that into a generalized histogram in the right panel. Errors in 
55(3839) were calculated as in Fig. |4] In this instance, we di- 
vide the CN-weak from CN-strong stars by shifting the base- 
line vertically until it encounters a significant gap. As in Fig. [4] 
CN-weak stars are shown as filled circles, and CN-strong stars 
as open circles. Although the two peaks in the right panel have 
roughly the same separation as in Fig. [4] the CN-weak group is 
a much more dominant component of the overall population. We 
have magnified the CN-strong curve in the right panel by a fac- 
tor of 30 (noted in the upper right corner of the panel) so that it 
is more easily comparable to the CN-weak curve. 

It must be noted that the slope of the baseline in Fig. 
is different from that in Fig. |4] This is unexpected, since (at 
fixed metallicity) the same processes ought to be moderating the 
progressive increase in CN bandstrength with increasing giant- 
branch luminosity: a declining temperature that both shifts flux 
from the science band of 5(3839) into the continuum band and 
permits more CN molecule formation, along with the "canonical 
extra mixing" decribed in Denissenkov & VandenBerg (2003) 
that progressively depletes carbon and enhances nitrogen abun- 
dance in the photosphere. We attribute the difference in base- 
line slopes to the different mass distributions of the two samples. 
Specifically, the M3 giants are an old, single-age population, 
with masses all around O.8M , while the field giants have the 
possibility of being considerably younger, and therefore more 
massive. A younger and more massive red giant, at fixed lumi- 
nosity, will undergo more rapid evolution, leaving less time for 
deep mix ing and surface abundance changes. Indeed, as is dis- 
cussed in Gilroy (1989), stars with masses greater than 2.2M 
evolve along the giant branch too quickly to even begin deep 
mixing. As a result, younger giants in the field will have weaker 
CN bands at a fixed M r , with the effect of pulling down the over- 
all baseline slope. 

To extend this type of analysis to the full final data set, we 
must carefully isolate [C/Fe] and [N/Fe] variations from the un- 
derlying metallicity, which strongly affects the appearance of 
double-metallic absorption features like /13883 CN. To that end, 
we divide the 1958 stars in the final data set into 0.1-dex-wide 
bins in [Fe/H] and search the bins independently for a bimodal 
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Fig. 7. CN versus CH bandstrength for the -1.4 < [Fe/H] < 
—1.3 metallicity bin, divided into three luminosity groups. Stars 
with both strong CN and relatively weak CH are shown as open 
circles. 



Fig. 6. Raw 5(3839) vs M r data for the eight 0.1-dex-wide 
metallicity bins in our sample. The baselines against which 
(55(3839) is measured are shown as dotted lines, and the max- 
imum metallicity for each panel is given in the upper left corner. 



distribution of CN bandstrength at fixed luminosity. Since CN 
bandstrength variations become very small at low metallicity, 
even for significant variations in [C/Fe] and [N/Fe] (e.g., Martell 
et al. 2008a; Briley et al. 1993), we limit our sample to the rela- 
tively metal-rich end of the halo. Figure|6]shows the raw 5 (3839) 
vs M r distribution for field stars in each of the metallicity bins 
between [Fe/H]= -1.0 and [Fe/H]= -1.8. The maximum metal- 
licity of the bin is given in the upper left corner of each panel. 
As in the left panel of Fig. [4] the dashed line in each panel is the 
baseline against which (55(3839) is measured. 

There are several features of note in this figure: the slope 
of the baseline flattens monotonically with dropping metallic- 
ity, from -0.11 in the -1.1 < [Fe/H] < -1.0 bin, to -0.04 in 
the —1.8 < [Fe/H] < —1.7 bin. I t is known that deep mixing is 
more efficient at low metallicity: iMartell et alj d2008d) find that 
drC/Fel/dMy is twice as la rge at [Fe/H]= -2.3 as at -1.0, and 
Sweigart & Mengel d 19791) predict less-compressed hydrogen- 
burning shells, allowing for more penetration by meridional cir- 
culation currents, at lower [Fe/H]. In addition, the gap between 
the CN-strong and CN-normal groups shrinks as metallicity de- 
clines, an effect that occurs even without reducing the size of 
variations in [C/Fe] and [N/Fe]. 

In order to identify stars with globular cluster-like carbon 
and nitrogen abundances in Fig. |6j we look for the classic 
CN-CH anticorrelation, used in globular cluster studies as a 
clear signal of carbon depletion and nitrogen enrichment. Figure 
Q sh ows d>5(3839) v ersus the CH bandstrength index 5 (CH) 
(iMartell et"afl l2008al) for the subset of the final data set with 
-1.4 <[Fe/H]< -1.3. The stars are further subdivided by abso- 
lute M r magnitude, since CH bandstrength is a sensitive function 
of both temperature and carbon abundance, and both decrease 
with rising luminosity. Stars with relatively large 55(3839) and 
relatively weak 5 (CH) in each luminosity sub-bin are shown as 
open circles. 

We repeated this selection process, dividing the metallicity 
bins into luminosity sub-bins, for the other seven metallicity 
bins, and Fig. [8] shows (55(3839) versus S{CH) for each of the 
metallicity bins, with all luminosity sub-bins collapsed together. 
Stars with strong CN (and weak CH, in the lower-metallicity 
bins), relative to the majority of halo field stars, are shown as 
open circles in each panel. Altogether, we identify 49 stars (also 



shown as open circles in Fig. [8]) as relatively CN-strong and CH- 
weak. Our CH bandstrength index 5 (CH) is calibrated based on 
bright red giants in the low-metallicity globular cluster M53. It 
is therefore less responsive to variations in carbon abundance in 
high-metallicity stars and in fainter giants than it was designed 
for. However, it was the most responsive over the full parameter 
range of our data, of the nine G-band indices we considered. As 
a result, we are less stringent in our CH bandstrength selection 
in the higher-metallicity bins than in the lower-metallicity bins. 

The eight panels of Fig. [9] correspond to the panels in Fig. 
[6] and show generalized histograms of (55(3839) for the CN- 
normal and CN-strong stars in each metallicity bin, using the CN 
bandstrength classifications made in Fig. [8] As in Fig. [5] the CN- 
strong curves needed to be amplified to be easily visible next to 
the CN-weak curves; the multiplicative factor is given in the up- 
per right corner of each panel, and the maximum metallicity in 
each bin is given in the upper left corner of each panel. These 
generalized histograms are qualitatively similar to Fig. [4] al- 
though with a clearly different CN-strong/CN-weak ratio: there 
is a separation between the two peaks of roughly 0.2 magnitudes 
in (55(3839), which shrinks as the overall metallicity drops. 

The presence of stars in the halo field with relatively strong 
/13883 CN absorption and relatively weak absorption in the 
/14320 CH G band, at fixed metallicity and luminosity, is a strong 
indication that globular clusters have contributed stars to the halo 
field. Given current models for the origin of light-element abun- 
dance variations in globular clusters, it does not seem possible 
that these stars formed in the halo with these atypical [C/Fe] and 
[N/Fe] abundances. An investigation of the abundances of O, Na, 
Mg and Al, which are also known to vary in CN-strong globu- 
lar cluster stars (see, e.g., Kraft 1994), would allow a more strict 
test, and possibly a more firm confirmation, of our claim that the 
49 halo field stars we identify here as CN-strong originated in 
globular clusters. 

For completeness, we also mention the possibility that these 
stars did not originate within globular clusters, but rather that 
their unusual abundances are the result of mass transfer from 
an AGB companion in a binary system. However, as is demon- 
strated in Lucatello et al. (2005), which discusses binarity in 
CEMP stars, the fraction of field stars expected to be in a bi- 
nary system with a companion of 3 - 8M Q and orbital parameters 
that permit evolution of the companion up to the AGB phase, and 
then mass transfer but not coalescence, is quite small. There is no 
high-precision radial velocity monitoring program in progress 
for these stars; we predict that such a program would be unlikely 
to find binary companions. 
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Fig. 8. CN versus CH bandstrength for all eight metallicity bins, 
selected from luminosity sub-bins as in Fig. [7] Candidate CN- 
strong stars are shown as open circles. 



4. Discussion 

Although light-element abundance variations have not been ob- 
served before in the halo, our identification of these CN-strong 
halo field stars is not wholly unexpected. There are several well- 
understood mechanisms for globular cluster mass loss, and the- 
oretical studies of globular cluster formation and evolution (e.g., 
D'Ercole et al. 2008; Baumgardt et al. 2008) predict significant 
mass loss in individual clusters as well as a dramatic reshaping 
of the cluster mass function with time. The data set we analyzed, 
selected from the SEGUE survey, is not representative of the 
full halo, in mass, evolutionary phase, or metallicity. Distances 
to the candidate CN-strong stars range from 4 to nearly 40 kpc, 
with 93% found within 20 kpc of the Sun. However, our result 
for red giants is generalizable to all halo stars, since abundance 
bimodality is observed to exist at all masses and evolutionary 
phases in globular clusters. More fundamentally, all of the cuts 
we made in selecting the final data set are blind to CN and CH 
bandstrengths and light-element abundances, and all nitrogen- 
enhanced giants with metallicities above [Fe/H] = -1.8 ought to 
show clearly strong CN bands. Since approximately 2.5% of our 
halo red giants exhibit strong CN bands and weak CH bands, we 
expect that the same fraction of the entire halo will contain the 
same abundance enhancements and depletions. This prediction 
can be confirmed by observations of dwarfs in the halo field, be- 
cause main sequence stars in globular clusters show the same 
abundance division as giants (e.g., Briley et al. 2002). 

In order to convert k, the present-day fraction of CN-strong 
halo stars, into f^ c , the fraction of globular cluster-originating 
stars in the halo field, we must consider what fraction CN-strong 
stars comprise of the stars originally formed in globular clusters. 
In the two-generation scenario of D'Ercole et al. (2008), roughly 
90% of stars originally formed in a globular cluster, consisting 
entirely of first-generation stars (with halo-like chemistry), are 
lost between the epoch of cluster formation and the present day. 
This means that /", the fraction of stars that remain as members 
of the globular cluster they were formed in, is around 0.1. Since 
/g P c , the fraction of present-day globular clusters stars that are 
CN-strong, is around 0.5 (e.g., Kraft 1994), we can calculate 



Fig. 9. Generalized histograms of 6S (3839) for CN-weak and 
CN-strong stars in each [Fe/H] bin. As in Fig. [6] maximum 
metallicity for each panel is given in the upper left corner. As 
in Fig. [5] CN-strong histograms were amplified for clarity. The 
multiplicative factor for each panel is given in the upper right 
corner. 



h -. Since V - 0.025 in the present study, this 



f^xfl ■'h 

t. _ _O025_ _ o 5 ; and that a rem arkable 50% of 



that if = 

suggests that /f = ^^ 
the halo field originally formed in the massive star clusters that 
were progenitors of the present-day globular cluster population, 
with a further unknown contribution of CN-weak stars made by 
globular clusters that did not survive to the present day and were 
not massive enough to self-enrich. 

While some numerical studies of galaxy formation (e.g., 
Boley et al. 2009) have suggested that the halo could be con- 
structed entirely from disrupted globular clusters, there is not 
presently a strong consensus on the role of globular clusters in 
cosmological-scale galaxy formation. Precise numerical study 
of the dynamical evolution of globular clusters is very com- 
plicated: the number of particles is large enough, and the rele- 
vant timescales short enough, that highly accurate simulations 
are very time-consuming. However, the development of semi- 
analytic prescriptions for the mass evolution of globular clus- 
ters would allow single-halo-scale simulations like those of 
Johnston et al. (2008) to include them as a source of halo stars, 
and to predict what fraction of the halo field ought to originate 
in globular clusters. 
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